t -^ cV via SUSY FCNC couplings in the 

unconstrained MSSM 

o 

S ! Jian Jun Liu°, Chong Sheng Li"* Li Lin Yang" and Li Gang Jin^ 

^ ' "• Department of Physics, Peking University, Beijing 100871, China 



p. 



b 



Institute of Theoretical Physics, Acadeniia Sinica, P. O. Box 2735, Beijing 100080, China 



m ; February 2, 2008 

in 
in 



Abstract 



o 
o. 

(— I ■ We recalculate the branching ratios for t —>■ cV {V = g, 7, Z) induced by SUSY 



FCNC couplings within the general unconstrained MSSM framework using mass 

Oh! 

O . eigenstate approach. Our results show that the branching ratios for these decays 

are larger than ones reported in previous literatures in the MSSM with R-parity con- 



^ ■ servation, and they can reach ~ 10 ^, ~ 10 ^, and ~ 10 ^, respectively, for favor- 

C^ I able parameter values allowed by current precise experiments. Thus, the branching 

ratios for t ^ eg and t ^ cj may be measurable at the LHC. 
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1 Introduction 

The top quark flavor changing neutral current (FCNC) processes t -^ cV {V = g,'y,Z) 
have tiny branching ratios in the standard model (SM) [1], and are too small to be mea- 
surable in the future colliders, and thus any detected signal of these rare decay events 
definitely indicates some new physics beyond the SM. Actually, t ^> cV {V = g, 7, Z) have 
been studied in various new physics models beyond the SM in detail, such as the two- 
Higgs-doublet model (2HDM) [1, 2], the technicolor model (TC) [3], the top-color- assisted 
technicolor model (TC2) [4], the models with extra vector-like quark singlets[5], the min- 
imal supersymmetry (SUSY) extension of the SM (MSSM) with R-parity conservation 
[6]- [11] and without R-parity conservation [12]. The decay branching ratios for t -^ cV 
{V = g,'~f,Z) are enhanced in general several orders of magnitude in these new physics 
models. The MSSM, which is believed as one of the most attractive candidates of new 
physics model, has gotten many attentions, and the investigation oit^cV{V = g, 7, Z) 
in the MSSM is a long story. C.S. Li et al. studied one-loop SUSY-QCD and SUSY-EW 
contributions in Ref. [6], subsequently G. Couture et al. recalculated and generalized the 
SUSY-QCD corrections to include the left-hand (LH) squark mixing in Ref. [7] and the 
right-hand (RH) squark mixing in Ref. [8]. All works above are within the framework 
of the MSSM with flavor-universal soft SUSY breaking terms. Later J.L. Lopez et al. 
further generalized the SUSY-EW corrections to the case of including neutralino-quark- 
squark loops in Ref. [9]. and CM. de Divitiis et al. reinvestigated t -^ cV {V = g,'~f, Z) 
in the universal case as well as non-universal case in Ref. [10], and obtained different 
results from Ref. [6, 7] due to the calculation of the relevant SUSY mixing angles and 
diagrams involving a helicity flip in the gluino line, which was conflrmed by J. Guasch et 
al. in a RG-based framework for t ^ eg decay [11]. All the above results of the MSSM 
are summarized in Table 1 (for comparing, we also list the results of the SM and new 
physics models mentioned above), one can flnd that they are all below 5 x 10"^, which 
is the roughly estimated sensitivities for the measurements of top rare decay at the LHC 
with 100 fb~^ of integrated luminosity [11]. 

However, all the previous works are limited to some constrained MSSM, in which some 
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Table 1: t -^ cV {V = g,'~f,Z) branching ratios of previous calculations. 'CKMl' and 
'CKM2' refer to models with extra vector-like up-type quark singlets and down-type quark 
singlets, respectively [5], and 'RPV refers to SUSY models allowing R-parity violation. 

strong assumptions or additional parameters besides ones in the MSSM are introduced to 
describe the FCNC couplings, but no any strong theoretical reasons of them have been 
found so far. It is necessary to study the FCNC top quark decays in the unconstrained 
MSSM [13], where the assumptions about the soft SUSY breaking terms are relaxed and 
new sources of flavor violation are presented in the mass matrices of sfermions, and conse- 
quently, some large contributions to FCNC processes induced by SUSY FCNC couplings 
(neutralino-quark-squark coupling and gluino-quark-squark coupling) can be obtained. 
Since the contributions to the top FCNC decays mediated by the charged current inter- 
actions (from W^, H"^, G"^ and x"^) are invisibly small as shown in the previous works 
[6]- [11] and can not be enhanced in this framework, in this paper we will reinvestigate the 
t ^ cV {V = g, 7, Z) only via SUSY FCNC couplings in the unconstrained MSSM, and 
try to show what are the maximal branching ratios for t -^ cV {V = g, 7, Z) in the MSSM 
using SUSY parameters allowed by current data, and whether they can be detected at 
the LHC. 



2 The t ^ cV {V = g,-f,Z) process induced by SUSY 
FCNC 

In the super-CKM basis [13], in which the mass matrices of the quark fields are diagonal 
by rotating the superfields, the up squark mass matrix 7W^ is a 6 x 6 matrix, which has 
the form: 

^ {Ml)^^ + {mlcos2(3m%[l-lsmHw))h (m|)^^ - /x(m„cot /^) II3 ^ 

^ {Ml) [^ - ^(m„ cot (3) h {Ml) ^^ + (ml cos 2/5 m| (i - | sin^ 9w) ) II3 j 

(1) 
where 9w is the Weinberg angle, II3 stands for the 3x3 unit matrix, the angle P is defined 

by tan/3 = f2/'i^i, the ratio of vacuum expectation values of the two Higgs doublets, fj, 

is the Higgs mixing parameter in the superpotential, and (M?.)^,^, {M^)jifi, and {M^)lji 

contain the fiavor-changing entries, which are given by 

respectively. Here Mq, M^^ and Au^, are the soft broken SU{2) doublet squark mass 
squared matrix, the SU{2) singlet squark mass squared matrix and the trilinear coupling 
matrix, respectively. They are directly related to the mechanism of SUSY breaking, and 
are in general not diagonal in the super-CKM basis. Furthermore, (M?)^/?, arising from 
the trilinear terms in the soft potential, namely Ajj^ijHjj UiUp is not hermitian. The 
matrix Ail can further be diagonalized by an additional 6x6 unitary matrix Zjj to give 
the up squark mass eigenvalues 

{MlY'''' = zlMpu. (3) 

Thus, we get new sources of fiavor-changing neutral current: neutralino-quark-squark 
coupling and gluino-quark-squark coupling, which arise from the off-diagonal elements of 
(M?)ll, {MI)lr and {MI)rr, and can be written as (/ = 1, 2, 3,i = 1, ..., 6, j = 1, 2, 3, 4) 

r - Qir - qis ■■ iV2gsT,^, [-{Zu)nPL + {Zu)(i+zy.PR 

—el 

X? - qir - qis ■ ^'^rs{[^= {Zu)lii-Sw{ZN)lj + Cw{Zn)2j) - F/(^t/)(/+3)j(^iv)4i]^L 

y ZswCw "^ 



+ [-^- (^l/)(/+3)j(^Af)li - Yu (^f/)/i(^Af)4j]-PR}- 

Here sw = sin6'vi/, cw = cos^^y, Pl,r = (1 Tln)/^, T"g is the SU{3) color matrix 
with color index a, r, s, and the unitary transformation Z^r diagonalizes mass matrix of 
gauginos and higgsinos to obtain the neutralinos. Thus the flavor changing effects of soft 
broken terms Mq, M^ and Ajj on the observables can be obtained through the matrix 

Zu. 

For the aim of this paper, the following strategy in the numerical calculations of the 
decay branching ratios of t ^ cV will be used: first we deal with the LL, LR, RL and 
RR blocks of the matrix Ai'^ separately and in each block we only consider the effects of 
individual element on the top quark rare decays, and then we investigate the interference 
effects between some different entries within one block and the interference effects between 
different blocks. In order to simplify the calculation we further assume that all diagonal 
entries in {M^)ll, {M^)lr, {M^)rl and (M? )/j/j are set to be equal to the common value 
MgugY, and then normalize the off-diagonal elements to M|usy [14, 15], 



yOu)LL- ^^ , [.du)RR- —rj2 



,ij^ _ V^^u)LL ..ij. _ V^^ij)RR 






i^u)LR = \7f^, {S}J)rl = ^tjP^, (^^j,^,j = 1,2,3). (4) 

'SUSY ^^-'SUSY 



Thus {M^)ll can be written as follows: 
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iMf,)LL = MlusY 



1 


{SF)ll 


iS}f)LL \ 


iSh')LL 


1 


K)ll 


mLL 


(sFUl 


1 / 



(5) 



and analogously for all the other blocks. 

The related Feynman diagrams for t -^ cV {V = g, 7, Z) induced by the SUSY FCNC 
are shown in Fig.l. Neglecting the charm quark mass, the amplitude of the decay process 
is given by 

M = u{p,)V^'u{pt)e,ik,X), (6) 



where Pt,Pc, and k are the momenta of the top-quark, charm-quark, and gauge boson 
respectively, and e^(A;, A) is the polarization vector for the gauge boson. The vertex V^ 
can be written as 

V^ = -^Y{PLF^, + PnF^,)-^^{PLF^, + PRF^,), (7) 

mt 

where FylJ) are the form factors, and their explicit expressions through the SUSY-QCD 
FCNC (^" -Qi- qi) are: 

(2Co'o + C'oiml - ml) - C^t " ^o) " SV^V^rBI] - Sm-gV.V^LVmB'o} , (8) 



g Qi 

si ~ -'si 



F^l = Filiy,L,R^V,R^L,VjL,R^V^R^L). (9) 



tJ- r-a 
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-m5\/7i?[(C^o' + Ct + C2')V^3 - HCl + C2^)\4i]}, (10) 

^2 = Fg',{V5L,R^V5R,L,VrL,R^VrR,L), (11) 

i^i = Y^^^J2i^t^'rL{'2C^oV3y5L + V;V,RB()+mgViVj^^^^ (12) 

F^^l = F^^i(\/5L,i?--V^5i?,L,VVL,i?--VVi?,L), (13) 

F^2 = F^2{V,L,R^V,r^l.Vjl,R^Vjr^l). (15) 



Fi2 = ]r^Y.^tViV,L[rnt{Ct, + Ci + C^,^V,L-rn-,{C^,+Ct + Cl)V,Rl (14) 



Fl, = Y^EKV^7L(E2C'o^oV^3'V^5L + KiV^5i^i?(: 



1=1 l'=l 



+mgVi'LV7R{V,LB'o - V^rBI;)], (16) 

^#1 = i^ll(V;'l--\/4'^,V^5L,iJ^V^5i?,L,V9L,i?--V9i?,L), (17) 

^#2 = 77^ E ^*V^3"^5L[m,(Ci^2 + C'2 + ^22)V^7L - ^^(^0^ + C^ + C^,)VrR], (18) 

^#2 = i^l2(V^5L,i?^V^5ii,L,V^7L,K^V^7i?,L), (19) 

and the explicit expressions through the SUSY-EW FCNC {xt ~ Qi ^ Qi) are: 

■rpa 6 4 

F,\ = T7-5-EE{2Co>*V^ilV^2lV3 + [(5o"^,o1^il 

lU/( lILt i^i i^^i 



+BlmtViR)V2R - B'om^oViRV2L]Vi}, (20) 

^,1 = F^\iV^L,R^V^R,L,V2L,R^V2R,L), (21) 

Fg2 = '^T.T.^tViLV^[{Ct2 + C^2+C'22)mtV2L-{C^o+Ct + C^)m^oV2R], (22) 
°^ 1=1 k=i 

F^2 = Fi2iyiL,R^ViR^L.V2L,R^V2R,L), (23) 

F^il = F;^{^{V,^Vi,V,-^ViT'^^l), (24) 

^zi = T77^JlJl{Jl'iCl^mtViLV2LV^ + V^L{-Blm^<.A^iRV2L + B\mtViRV^ 

iOTT Ult i^^ ^^^ i,^-^ 

4 

+B^m^^>^ViLV2R) + Y. mtViL[{mlV2LV8R - m^o{mtV2R + m^oy2L)V8L)C^ 

k'=l 
+ V2lV8rB^ - 2V2lVsrC;^o + {m]V2LV8R + mtm^(>y2RVsR 

-mtm^oV2RVsL)C:^]}, (25) 

Fz2 = ^i:T.mtViL{J:V^[{Ct2 + C^ + C^22)mtV2L-{CS + Ct + C^)m^oV2R] 

°^ 1=1 k=l V=l 

4 

+ E [-m^iy2RVsLCl - V8R{mtV2LCl2 + {mV2L + m^o^y2R)C^2 + 
fc'=l 

mtV2LC'22)]}^ (26) 

Hereof = 5,(0,m|o,m| ), B^ = 5,(m2,m|„,m? ),5f = fi,(0,m|o,m|o;, Bf = 5,(0, m?,m?), 
5<^ = 5i(0,m?,m|), 5/ = 5i(m2,m?,m|), C«.^. = ^,^^(0, 0,m2,m|o,m? ,m? ), 

and Cj^jj = Cj^jj(0, 0, m^, m|o , ""^l , ""^l , ) are 2 and 3-point one-loop integrals [16]. And the 
relevant couplings are: 

— e 1 

^IL = '^{'^^ {Zu)2l[-Sw{ZN)lk + Cw{ZN)2k]—yu{^u)bl{ZN)ik}-, 

yZSwCw "J 

2^26 

VlR = "^[-^ {Zu)5liZ]\r)ik — Y^{Zu)2l{Ziy)4k], 

oCw 

— e 1 

V2L = i{~/^ iZu)3i[-Sw{ZN)ik + Cw{Ziy)2k] —yu{Zu)Qi{Ziy)4^l:}, 

y/ZSwCw "J 

V2R = i[— {Zu)6l{ZN)ik — Y^{Zu)3l{ZM)4k], 

oCw 



ZSwCw /=! <J 

VsL = -iV2gs{Zu)2h V^R = i^/2gs{Zu)-,u Vq = ~igs, 

Vtl = -iV2gs{Zu)3i, Vjr = iV2gs{Zu)&u 
g 

y^L = i- [(^Af)4fc(2'Af)4fc' — (2'Af)3fc(2'Af)3fc/], 

ZSyi/Cw 
g 

VsR = —i- [(^7v)4fc(^7v)4fc' — (^7v)3fc(^7v)3fc']- 

After squaring the decay amplitude and multiplying by the phase space factor, one obtains 
the decay width oi t ^ cV {V = g, 7, Z): 

Tit ^ eg, C7) = ^m*[(2i^yi " ^yi)^yi* + (2i^yi - ^y2)i^yr - {F^i + F^2)F^; 
yoTT 

-iF^i + Fv2)F^2l (28) 

r(t - cZ) = ^Q^J^2^^M - mlf{2mlFll[ml{AFi, - F|,) + m?(2F|, + F|,)] 
+2m^F|*Kz(4F|i - F^) + m?(2Fi + F^)] - K - ^Dl^^i* 
(m|F|2 - m^(2Fi + F^)) + Ff^^lml^l^ - m2(2F|i + FI^))]}, (29) 

and we define the branching ratio as Ref. [1]: 

which will be the main object of our numerical study. 



3 Numerical calculation and discussion 

In our numerical calculations the SM parameters were taken to be rrit = 174.3 GeV, 
Mw = 80.423 GeV, Mz = 91.1876 GeV, sin^ ^vk = 0.23113 and ^.(M^) = 0.1172 [17]. 
The relevant SUSY parameters are /i, tan/5, Msusy and rrig, which are unrelated to 
flavor changing mechanism, and may be fixed from flavor conserving observables at the 
future colliders. And they are chosen as follows: Msusy = 400, 1000 GeV, tan/3 = 4,40, 
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rrig = 200, 300 GeV and fi = 200 GeV. As for the range of the flavor mixing parameters, 
{Su)ll are constrained by corresponding ((5£))ll[15, 14, 18, 19], in which {5\^)ll also is 
constrained by the chargino contributions to K-K mixing[20], and Dq-Dq mixing makes 
constraints on {5}f)LL, {5}f)LR and {5}f)RL[2l]. And {51^)ll, {5f)LL, iSf})RL and i6fj^)RL 
are constrained by the chargino contributions to B^-Ba mixing[18]. Finally, there also 
are constraints on the up squark mass matrix from the chargino contributions to b ^ 
S7[22, 14]. Taking into account above constraints, in our numerical calculations, we use 
the following limits: 

(i) {61j')ll, {Sh')LR and {6h')nL is less than 0.08Msusy/(1 TeV); 

(ii) {61j^)rr and {61j^)ll are limited below 0.2Msusy/(1 TeV); 

(iii) i6ff)LL, {5lf)LR, {5lf)RL, {5lf)RR. {5}f)LR, {^W) RL and {5}f)RR Vary from to 1. 

First of all we should point out that the contributions arising from SUSY-EW FCNC 
are in general at least one magnitude of order smaller than ones arising from SUSY-QCD 
FCNC, and the dominant contributions to the decay branching ratios come from the 
latter. Furthermore, our calculations show that the decay branching ratios only weakly 
depends on tan/?, so we only discuss the results in the case of tan/3 = 40 below. Our 
results are shown in Figs. 2-3, where there are three common features of these curves: the 
first is that the branching ratio increases rapidly with the mixing parameters increasing, 
and the second is that the branching ratio depends strongly on the gluino mass m^, and 
the last is that the dependence of branching ratio on the Msusy is medium comparing 
with above two parameters. 

For each decay modes t -^ cV {V = g, 7, Z), in Figs. 2 and 3 we show the dependence 
of the decay branching ratios on RR and LR off-diagonal elements, respectively (We do 
not show the results for LL off-diagonal elements as their contributions are similar to the 
ones for RR off-diagonal elements). We find that for t -^ cV the largest results come 
from the LR block, which arises from the soft trilinear couplings Au. We also give the 
results of the interference effects on the branching ratios for t -^ eg in Figs. 2(c)and 2(d). 
In general, these interference effects increase the decay branching ratios, and since the 
interference effects between different blocks are similar, we do not show them for the space 



of this paper. The results of decay t -^ cy, cZ are about two orders of magnitude smaller 
than ones of decay t -^ eg, as shown in Figs. 2 and 3. From these figures, we can find 
that the decay branching ratios for t ^> cV {V = g,'y,Z) induced by the SUSY FCNC 
couplings can reach ~ 10^^, ~ 10^^ and ~ 10^^, respectively, for the favorable parameter 
values allowed by current precise experiments, and they are larger than all the previous 
ones in the MSSM with R-parity conservation (it should be pointed out that the results 
of Ref. [6], [7] and [8] in Table 1 are obtained at rrig = 100 GeV, which is disfavored by 
current data). 

According to the analysis of T. Han et al. [23], the sensitivities for t —^ c-f and t ^ cZ 
at the LHC with 100/6~^ integrated luminosity are 5 x 10~^ and 2 x 10~^, respectively, 
and our results show that the rare decay t ^ cf may be detectable. Later T. Han et al. 
[24] and M. Hosch et al. [25] studied the sensitivities to the top quark anomalous FCNC 
couplings at the LHC for single top quark and direct top quark productions, respectively, 
and the corresponding decay t ^ eg branching ratios transferred from their results are 
4.9 X 10~^ and 2.7 x 10~^, respectively. Thus, our results of the branching ratios for 
t ^ eg indicate that the top quark FCNC production processes (both for single top and 
direct top) may be measurable at the LHC. But if we use the 5 x 10^^ as the sensitivity 
for the FCNC decay t -^ eV at the LHC with 100/6"^ integrated luminosity as shown in 
Ref. [11], our results show that the rare decay t ^ eg are also potentially measurable at 
the LHC. 

In conclusion, we have calculated the top quark rare decay t -^ eV {V = g,'~f,Z) 
induced by SUSY-FCNC couplings in the general unconstrained MSSM using mass eigen- 
state approach. Our results show that the branching ratios for these decays are larger 
than ones reported in previous literatures in the MSSM with R-parity conservation, and 
especially, the branching ratios for the rare decay modes t — * eg, 07 we calculated are 
very hopefully to be measurable at the LHC for the favorable parameter values allowed 
by current precise experiments. Moreover, we find that the decay branching ratios for 
t -^ eV(y = g, 7, Z) strongly depend on the soft trilinear couplings Au, and it is possible 
to get some valuable information of soft SUSY breaking parameters by measuring the 
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branching ratio for the top quark rare decay at the LHC 
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Figure 1: Feynmann diagrams for t -^ cViV = g, 7, Z). 
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Figure 2: The decay branching ratios for the t ^ eg with mixed RR off-diagonal el- 
enients(a) and LR off-diagonal elements(b), and the typical interference effects of RR 
block(c) and LR block(d). Here, sohd line: m-g = 200 GeV, MgusY = 400 GeV; dashed 



hne: m-g = 300 GeV, MgusY = 400 GeV; dotted line: m-g = 200 GeV, MgusY 
GeV; dash-dotted line: m^ = 300 GeV, MgusY = 1000 GeV. 



1000 



15 




o 
m 



o 
m 




Figure 3: The decay branching ratios for the t —* c-f with mixed RR off-diagonal ele- 
nients(a) and LR off-diagonal elenients(b), and the decay branching ratios for the t ^ cZ 
with mixed RR off-diagonal elements(c) and LR off-diagonal elements(d). Here, solid 



hne: m-g = 200 GeV, MgusY = 400 GeV; dashed line: m-g = 300 GeV, MgusY 



400 



GeV; dotted line: m-g = 200 GeV, MgusY = 1000 GeV; dash-dotted line: m-g = 300 GeV, 
MsusY = 1000 GeV. 
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